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Abstract 



O 

' We consider the eigenstate problem for a Hamiltonian operator of the field 

theory. Methods of construction the effective field theoretical Hamiltonians for 
, which the eigenstate problem may be solved are discussed. In particular, we discuss 

the method of flow equations from a general perspective as well as in application 
to the gauge field theories. Flow equations transform the Hamiltonian to a block- 
diagonal form with the number of particles conserved in each block and thus reduce 
the original bound state problem to a set of coupled eigenstate equations with an 
effective Hamiltonian in each sector. 

Applications of flow equations to the Hamiltonians of QED and QCD in the 
light-front gauge and the QCD Hamiltonian in the Coulomb gauge are considered. 
Using flow equations, we derive the effective Hamiltonians as well as the renor- 
CS| ■ malized gap equations and the Bethe-Salpeter equations for the bound states in 

I these theories. We show that the obtained equations are finite in both UV and IR 

' regions and are completely renormalized in UV, i.e. the corresponding solutions do 

not depend on the cut-off A. 



Qh. We calculate positronium spectrum, glueball masses, tt — p mass splitting, gluon 

^ ' and chiral quark condensates and compare our results with the covariant calcula- 

tions and experimental results. Use of flow equations to calculate the dynamical 
terms is critical to achive good agreement with experimental results. 



^Based on talks given at JLAB (August 30, 1999), Duke University (October, 1999), BNL (November 
3, 2000), MIT (December 11, 2001). 
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1 Flow equations. Idea and technique. 



1.1 Hamiltonian bound state problem 

Motivation 

Nonabelian gauge field theories, like Quantum Chromo dynamics (QCD), are well 
understood in the high-momentum region, since due to the asymptotic freedom cou- 
pling constant is small and Feynmann rules of the covariant perturbation theory, using 
Lqcd, provide convincing agreement with experiment. However, low energy QCD, where 
such nonperturbative phenomena as confinement and chiral symmetry breaking are tak- 
ing place, is still understood only on a phenomenological level. Except for the lattice 
studies and some phenomenological models, there is a lack of analytical nonperturba- 
tive methods. Here a systematic analytical method to solve bound state problem, using 
Hamiltonian, 

HqcdI^I^) = E\i;), (1) 

is presented. 
Problems 

Two general problems occure when solving the eigenstate equation for a relativistic 
Hamiltonian of the field theory: field theoretical Hamiltonian is an infinite dimensional 
matrix, first, in the particle number space and, second, in the energy space. 

First, the number of particles is not fixed in the field theory due to allowed creation 
and annihilation processes in vacuum. Therefore, any physical state contains, in principle, 
infinite many Fock components 

\lp)meson = Cgg\qq) + Cqqg\qqg) + . . . . (2) 

Thus, we need a method to construct an effective Hamiltonian, which acts in a smaller 
truncated space and provides the same eigenvalues as the original Hamiltonian. 

Second, there are states with (infinite) large energies/momenta in the Hamiltonian, 
which contribute UV divergent matrix elements. Hamiltonian matrix in a momentum 
space (plane wave space) reads 





H = 



(3) 

A ^ oo 



where the top left matrix element correspond to small (zero) momentum, and the bottom 
right - to (infinite) large momentum with the cut-off UV scale Auv- Thus, we need a 
method of UV regularization and subsequent renormalization directly for a Hamiltonian 
operator. 

1.2 Hamiltonian methods 

There are following Hamiltonian methods, which address both problems. However, there 
are some disadvantageous in these methods, restriting the range of their applicability. 
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Figure 1: Projection of high Fock components to the low ones. 

1.2.1 Tamm-Dancoff truncation and iterated resolvents 

The method of Tamm-Dancoff truncation, treating the 'infinities' in Fock space, has 
been invented in 50's, and later generalized in the method of iterated resolvents by Pauli 
in 80's. The main idea of this method is to project high Fock components onto the 
lower ones in sequence, ending up with an effective Hamiltonian which acts in a smaller 
truncated space (Fig. 0) and therefore it can be solved for eigenstates and eigenvalues. 

Technically, by projecting high Fock states one inverts a resolvent, which is the energy 
denominator containing matrix element of the original Hamiltonian, {E — {n\H\m))~^ . 
Inversion is done within some approximation scheme. It turns out that this scheme 
breaks down for QCD, convergence cannot be achieved, and eigenvalues of an effective 
QCD Hamiltonian do not predict QCD mass spectrum correctly. Physically, starting 
with a canonical QCD Hamiltonian and working in terms of bare fields, one fails to find 
a representation for the QCD Hamiltonian with fixed number of particles, since QCD 
vacuum is unstable and due to a strong QCD interaction it creates and annihilates bare 
particles. In addition, unregulated UV divergences, comming from high energies, appear 
in the effective Hamiltonian. A way out to proceed with QCD is suggested by various 
phenomenological studies, where instead of bare (current) degrees of freedom one uses 
the renormalized (constituent) ones. In other words, one needs a scheme to regulate and 
renormalize Hamiltonians. 

1.2.2 Similarity renormalization 

Several years ago, in 1994, the similarity renormalization scheme, which treats 'infini- 
ties' in the energy space, has been suggested by Glazek and Wilson. Similarity scheme 
is an analog to the Wilson's renormalization through an effective action where high en- 
ergy modes integrated out in a path integral, but formulated directly for a Hamiltonian 
matrix. The aim is to regulate and renormalize away UV divergencies. The main idea 
of the similarity renormalization is to find a unitary transformation, U, which brings a 
Hamiltonian matrix to a band-diagonal form (Fig. |^) with the width of the band. A, 
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Figure 2: Band-diagonal Hamiltonian in the energy space. 



being much smaller than the UV cut-off, Aj 



U\X,A)H{A)U{X,A) 
\Ei - Ej\l\ < A . 



^e//(A,A) 



(4) 



High and low energy states decouple from each other in a band-diagonal effective 
Hamiltonian. Indeed, high order in the perturbation theory 



with A <C A, connects high and low energies. Therefore, a low-energy effective Hamilto- 
nian, i^e//) can be taken and solved for few lowest states. Since transformation is unitary, 
eigenvalues of i?e// should be the same as low lying states of the original Hamiltonian. 
Similarity renormalization is formulated as a perturbation expansion in coupling con- 
stant. As one shrinks the band width, A, one breaks perturbation theory at some point, 
since diagonalization is a nonperturbative problem. In similarity, higher orders of pertur- 
bation theory correspond to higher Fock states. Therefore, breaking of the perturbation 
theory means that mixing of higher Fock components becomes important and cannot 
be neglected. In QED, a perturbative mixing of electron-positron and clcctron-positron- 
photon states provides reliable approximation for positronium bound states. Also, the 
window for the band width can be easily found, ma^ -C A -C ma where a = e^/47r. 
However, for QCD a perturbative mixing does not work. One needs a nonperturbative 
scheme to incorporate high Fock states mixing. 

1.3 Flow equations 

Method of flow equations for Hamiltonians was suggested independently by Wegner in 
1994. Apphed to the field theory, flow equations try to solve both problems: 'infini- 
ties' in the Fock space as well as in the energy space. H^jf incorporates effects from 
high Fock components and large energies. This method can be viewed as a synthesis 
of the Tamm-Dancoff and similarity renormalization approaches. However, instead of 
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,2A/A ^ ^ ^ 



(5) 
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Figure 3: Block-diagonal Hamiltonian in the particle number space. 



band-diagonalization in the energy space, that is hard to achieve, one block-diagonalizes 
Hamiltonian in the particle number space (Fig. |^). 

One seeks for a unitary transformation, which brings a Hamiltonian matrix 

^ (PHP PHQ\ 

[ QHP QHQ J ' ^''^ 

where P and Q = 1 — P are projection operators in the particle number space, to a 
block-diagonal form 

" ^ ( QH,,Q ) • (^) 

with the number of particles being conserved in each block. Since blocks decouple from 
each other in a block-diagonal effective Hamiltonian, full bound state problem is reduced 
to several eigenstate equations in each block. Usually it is easier to diagonalize separately 
P or Q sectors than to diagonalize the original Hamiltonian matrix. 

Unitary transformation is governed by a continuous flow parameter /, 

U{l)HU-\l) = H{1), (8) 

with initial condition [/(/ = 0) = 1. As / — > cxd an effective Hamiltonian coincides with 
its diagonal part, Eq. (|^), H{1 oo) = PH^ffP. Unitary transformation is given by 

U{l)=Tiexp(^j\{l')dl'^ , (9) 

where ri{l) is the generator of transformation, and Ti is the /-ordering. Wegner's flow 
equations are written in a diffrential form 

^ = m),Hii)] 

r/(/) = [H,{l),H{l)], (10) 
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where the generator, chosen by Wegner as a commutator of a diagonal part 

QH{l)Q 



m) - ( ] . (11) 



and an off-diagonal part, 

mi) - HAD = ( QHiDP ^'"'^ ) ■ (12) 

insures that an effective Hamiltonian is diagonal at I — > cxo, H{1 —>■ oo) = H^ll — > oo). 
Diagonal - particle number conserving, Eq. (|11|), and off-diagonal - particle number 
changing, Eq. ([T^) , parts of a Hamiltonian contain sector Hamiltonians which are fuctions 
of the flow parameter /. 

Using Eqs. (|10D and (|ll|)7 dH); flow equations for the diagonal and off-diagonal Hamil- 
tonians, and the generator of transformation are written in the operator form as 

^PHP = PriQHP - PHQriP 
dl 

^PHQ = Pr]QHQ- PHPr]Q 
dl 

PriQ = PHPHQ- PHQHQ . (13) 

Choosing eigenstates of initial P and Q sector Hamiltonians for a basis (suppose it can 
be done), Eq. (O) is given in the matrix form 



VpS) = {E,{l)-ES))hM, (14) 

with Pif(/)P Ep{l), QH{l)Q Eq{l), a.nd PH{l)Q hpg{l);p,p' are energy /momentum 
indices in the P space. Substituting the generator of transformation, we get a coupled 
system of equations for particle number conserving 

dhppi{l) _ / dhpq{l) hqpijl) hpg{l) dhqp>{l) \ 

dl ~ q\ dl Ep{l) - Eg{l) Ep,{l) - Eg{l) dl J 

and particle number changing sectors 
dhpg{l) _ 



dl 



{Ep{l) - Eg{l))%g{l) , (16) 
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respectively. Formal solution of these equations in P space is given at the flow parameter 
I ^ oo. Particle number conserving part 



VM = V'(0)- dlY. 7\ r'Ln. + 



^ ^ dl Ep{l)-Eg{l) Ep.{l)-E,{l) dl 



EAD - i^.(o)-///T ^^(,)!^^(,) |(MOMO), (17) 

includes the initial P sector Hamiltonian plus terms, generated by flow equations when 
eliminating particle number changing sector. If p corresponds to one-particle degree of 
freedom, the second equation in ([T7|) is the energy gap equation, and the first equation 
is the Bethe-Salpeter integral equation for a bound state. Eq. ( [T7| ) resembles the second 
order perturbation theory, where hpq{l) is a vertex and l/{Ep{l) — Eq{l)) is the energy 
denominator. It is important, that energies also fiow with /, and are given by diagonal 
matrix elements in P space, p = p'. To find P space effective Hamiltonian, hppi, one 
should know / dependent solutions for the energies, Ep{l) and Eg{l), and for the particle 
number changing part, hpq{l). Solution for the particle number changing part 



hpq 



(/) = hpq{0)exp(-jyi'{Ep{l')-Eq{l'))A, (18) 



characterizes by an exponential function. For /-independent energies, hpq decays expo- 
nentially to zero with I oo. Renormalization of energies, i.e. /-dependence of energies, 
brings some additional factor before the exponential decay. Energy renormalization is 
crucial for degenerate matrix elements, when initially \Ep{0) — Eq{0)\ = for / = 0. In 
this case, /-dependence of energies provides power law decay of hpq instead of exponen- 
tial one. Thus, the Wegner's generator (Eq. (^)) insures that the off-diagonal part is 
eliminated always, even for degenerate states. Eqs. (0) and (p^ ) are coupled, therefore 
they should be solved selfconsistently. 

Our strategy is to solve first the equation for an effective renormalized energy (the 
gap equation), and then using energy solutions to find an effective Hamiltonian in P- 
sector, hppf, (kernel of the Bethe-Salpeter equation), which is diagonalized numerically 
for bound states. Success of the procedure depends on a particular choice of diagonal 
and off-diagonal Hamiltonian parts, that is discussed in application to QED and QCD 
further. 

In a block-diagonal Hamiltonian, there is some freedom left to transform inside of 
each block without changing a block-diagonal structure. This freedom can be used by 
choosing different similarity functions, f{x), 

hpqil) = hpq{0)f{Xpq) 

Xpqil) = f dl\Ep{l')-Eq{l'))\ (19) 

with general properties to be equal unity at the origin, and to fall down at large arguments 
/(O) = l,/(x-^oo) = 0. (20) 
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particle number 



energy 



Figure 4: Pentadiagonal form of Hamiltonian in the particle number space; each sector 
contains matrix elements with all possible energy changes. 

Similarity function reflects the rate how quickly the particle number changing sector is 
eliminated. The generator is written through the similarity function as 

fpq{x) = f{xpg). Flow equation scheme with P and Q sectors can be easily generalized 
to the case of many Fock components (Fig. ^). It is important, that the relevant degrees 
of freedom are renormalized during this procedure. As a consequence, it is not coupling 
constant which is used as a small parameter, but rather some dimensionless combination, 
say Xpg = \hpg\/\Ep - Eg\. 

Elimination of the particle number changing sectors is carried out not in one step but 
sequentially for energy differences which exceed the band width A (flowing cut-off), Eq. 
®, 

A = ^ < \Ep{0) - Eg{0)\ < -^=L= = A ^ oo , (22) 

for /-independent energies, with / = (Fig- U)- 

Lowering the cut-off A, we effectively scale the theory down to low energies, that 
provides a connection between flow equations and similarity renormalization. System of 
coupled flow equations contain equations of the renormalization group, in addition to 
appearing equations for noncanonical operators. In this way, we incorporate renormal- 
ization in the many-body technique of block-diagonalizing. 

Summary I 

1) Summarizing, flow equations perform Hamiltonian 'renormalization' in the energy 
and particle number space, in the sence that an effective Hamiltonian is flnite in both 
spaces and contains effects of high energies and high Fock components. Technically, flow 
equations allow to include directly the dynamical perturbative corrections, responsible 
for the renormalization and the UV asymptotic region, into the many-body calculations. 

2) For the flrst time, due to the dynamical interactions, the physical equations are 
obtained UV and IR flnite. 
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Figure 5: Flow equations block-diagonalize the bare regulated by A Hamiltonian {H{A)) 
in the particle number space. For the finite A (after the unitary transformation U{A ^ 
oo, A) is performed) the matrix elements of the particle number changing sectors are 
squeezed in the energy band \Ei — Ej\ < A on the left hand side the picture and are 
eliminated completely as A — ^ (that corresponds to U{A — > oo, A 0)) on the right 
hand side of the picture. Block-diagonal effective Hamiltonian is also UV renormalized 
by flow equations. 
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2 Application of flow equations to the gauge fleld 
theories 



Flow equations have been apphed mostly in solid state physics by Wegner and his group. 
Also some toy models have been considered in the field theory. Here we consider several 
examples in the gauge field theories, such as the light-front QED [|l| and light-front QCD 
Hamiltonians as well as Hamiltonian of gluodynamics [§] and QCD Hamiltonian with 
dynamical quarks in the Coulomb gauge. In Hamiltonian dynamics of the gauge field 
theory, first a gauge should be fixed and constraint equations should be solved. Hamilto- 
nian operator is expressed only through the physical degrees of freedom. Generally fewer 
equations of motion and current algebra relations (like Ward identities) occure in the 
Hamiltonian dynamics than for an effective action. Moreover, there are no ambiguities 
in these equations since the physical fields have been used. However, one should keep 
track on covariance and gauge invariance manifest for physical observables. This serves 
to check the validity of the calculations and approximations. 



2.1 Flow equations in the light-front QED and QCD 

Light-front dynamics is similar to the equal time one, except that the light-front time is 
given by = t + z, and the commutation relations and initial conditions are formulated 
at a quantizaton plane = 0, instead of t = in the equal time framework. Evolution 
(propagation) is considered along the light-front, < x"*". In the collinear limit (small 
light-front x) which corresponds to the light-front origin point, the light-front field theory 
has a singular behavior associated with nontrivial vacuum effects in the equal time. 
The Lagrangian density for QED/QCD 

£ = ~F^,F^'' + i^{t^ + g4-m)i; (23) 

is considered here in the light-cone gauge A'^ = A^ + A'^ = 0. Zero modes will be 
disregarded. The constrained degrees of freedom, A~ and (A± = |7°7^ are projection 
operators, thus ip± = A±il), and ip = ip^ + ipJ) are removed explicitly and produce the 
canonical QED/QCD Hamiltonian. It is defined through the independent physical fields 
A±_ and '0+- To solve the constrained equations for A~ and ifj- the auxiliary fields 

A - A - ^ T+ 

^ = ^+ + (m/3 - ta'd^i) , (24) 

are introduced. The fermion current is J'"(x) = ^7^^. The resulting canonical Hamilto- 
nian = P+ is given as a sum of the free Hamiltonian and the interaction 

H = P+ = Ho + V + W . (25) 

The free Hamiltonian Hq, the quark and gluon kinetic energies, is 

Ho = lJ 7-^^^^^^Jj^^ + i^(.V^)^i,) . (26) 
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In the interaction V + W, the vertex interaction V is the hght-cone analogue of the 
minimal coupling interaction in covariant QED/QCD and W = Wi + W2 is the sum 
of the instantaneous-gluon Wi and the instantaneous-quark interactions W2- The latter 
arise from the constraint equations (analog of the hght-front Gauss law). More exphcitly, 
the interaction is given by 

2 J [io+) 



2 + 

W2 ^ ^1 dx+dw ^rA^ (7^^*) , 



v 



I dx+dWB^'B;^. , (27) 
where the current includes the quark and gluon parts 

>(x) = ¥7'**^"+-^F'^^^], (28) 

with — A^T"-. The instantaneous gluon and quark interactions behave as and 
l/g"*", respectively, with the momentum transfer q — {q±, g+). 

It is convenient to work in a second quantized form, decomposing the physical fields 
through creation and annihilation operators. By definition, the fields ^ = ^+ + and 
A'' = (0,lx,i+) are the free solutions which in momentum space are parametrized as 

A ^2p+(27r)3 ^ ^ 

M^) = T. I (a(p)6,(p,A)e-^- + at(p)6-(p,A)e+-^-) . (29) 

A •' y2p+(27r)"^ 

The single particle operators obey the commutation relations 

'a{p),a\p')] = {b{p),b^{p')} = {d{p),d\p')} = S{p+-p+')S('\p^-pi)S^' . (30) 
The free part is given by 

Ho{l) = / ^ E{p; I) {h\p)h{p) + S{p)d{p)) + / u;(p; a\p)a{p) .(31) 

The single particle energies {E = p~) depend on the 3-momentum p — {p'^,p±) 

pl + m\p]l) P±+f^\p]l) 
-^(^' ^ ^+ ' ^^P' ^' ^+ ' (^^^ 

and potentially on the flow parameter through the mass m'^{p] I) of the particle in ques- 
tion. Inserting the free fields into the Hamiltonian yields for the quark-gluon vertex 
interaction 

m ^ '[\^[\^ 5i^\p.-P2-P2) 



7(2^^ ^/2pt■^ ^J^J ^2^ 
[g{Pi,P2,Pi\l)i>^{pi)Kp2)a{pz)u{pi)i{pi)u{p2) + ..^ , (33) 
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and the quark- antiquark interaction in the exchange channel 



1 



(27r)3 
W, 



[d^Pi] [[d'p2] [[d^Pz] f[d^P4\ 



2pi 



5^^\pi -P2+P3 -Pi) 



,iu{PuP2,P3,P4; I) b^{Pi)b{P2)d{-p3)d^{-p4)T^^T^4^ 

U{pi)j^u{p2)v{-P3)j''v{-P4) + ...] . 

The integration symbols denote 
[d^p] = dp+d^p^ ' 

X 



(34) 



(35) 



and the abbreviations u{p) = u{p, X) and = jf^efi^^p, X) are introduced. When uni- 
tary transformation with the finite fiow parameter / is applied, all coupling constants 
and masses become unknown functions of momenta and I (This is refiected in explicit 
dependence of g and m, /i, and Wfj,,^). Initial conditions are given at / = 0. The effective 
coupling constant has the initial value 



9(1 = 0) 



9- 



(36) 



with the fine structure constant a = (?^/47r ~ 1/137 in QED, and the strong coupling 
g = Qs iTi QCD. The two-point interaction includes the instantaneous and dynamical 
generated by fiow equations interactions, W = W^^'^^ + W^"^", with initial conditions 

W^-*(/ = 0) = W^-(/ = 0) = 0, (37) 



where the momentum transfer is q = pi — p2- 



2.1.1 Effective Hamiltonian for the light-front QED 

Representing through creation and annihilation operators, we define explicitly the par- 
ticle number conserving and particle number changing parts which directly define the 
generator of transformation. In QED, all canonical interactions are shown in the Table 
for example, the electron-photon minimal coupling is the matrix element H12, the electron 
and photon instantaneous interactions are ifi4 and H22, respectively. Particle number 
changing terms are off-diagonal in this table and are eliminated using flow equations. 
The generator of the unitary transformation is 

fj{i) = m)+ui)- (38) 

where the generator fji = [Hq, V] eliminates the electron-photon coupling V and is given 
by 

m-^fK^ fK^ fK^s^^\p.-P2-P3) 



(27r)3 J Pp+ J Pp+ J Pp\ 

'niPuP2,P3;l)blb2a3u{pi)^{p3)u{p2) + ... , (39) 
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Table 1: QED canonical Hamiltonian in the light-front gauge. 
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Table 2: Effective QED Hamiltonian generated through the second order by ffow 
tions in the hght-front gauge. 
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and the generator fj2 = [Hq, W] eliminates the off-diagonal two-point interactions and 
its explicit form depends on the particular sector where is W. 

QED coupling is a small parameter. Eliminating the particle number changing sector 
to the second order in coupling, we generate new terms which are given in the Table ^. In 
the second order 0{g'^), new dynamical two-point interactions arise from elimination the 
electron-photon vertex terms, while elimination of the instantaneous terms W generates 
new interactions in the third and high orders which are not depicted. In the Table ^ 
black area in the two-point interaction depicts an effective kernel which depends on all 
four in- and out-going momenta. 

Effective electron-positron interaction (matrix element H22) includes the instanta- 
neous and dynamically generated by flow equations terms. In the exchange channel, it 
is given by a product of the current-current term and the interaction kernel. 



where the photon interaction kernel 



B 



flU 



(40) 



(41) 



contains the integrals Ji and I2 which are defined by two similarity form factors in each 
vertex 



/i 



^^^Ql^)^^^^'' ^ 



(42) 



and for I2 the indeces 1 and 2 are interchanged. Here, and Q2 are momenta transfers 
in two verteces. In the light-front framework, they are given by 



Ql 



{x'k± — xk'j_) + m (x — x') 



/\2 



XX' 



Ql 



',2 — Ql\x-*{l-x); x'^{l-x') ■, (43) 

On mass shell, these momenta are reduced to the three photon momentum transfer, 
-12 /n2 , zfi 



(44) 



Using three choices of similarity function; exponential / = exp (— Q^/A^), gaussian / = 
exp(— Q^/A^) and sharp / = ^(1 — Q^/A^); we obtain the following kernels for the 
effective electron-positron interaction. 



B 



B 



B 



1 

1 
1 



9iJ.v ^2 



9iJ.u ^2 



9tJ.v ^2 



9^lu 






q+2 


9f,u 






q+2 



5Q' 



(45) 
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respectively. Here, is the average momentum transfer and dQ"^ shows the off-shellness. 
On the energy shell, i.e. SQ"^ = 0, the electron-positron interaction is reduced in each 
case to the 3-d Coloumb potential with behavior l/(f. 

Picking the ee-sector, we solve an effective interaction for positronium bound states. 
The light front Schrodinger equation for the positronium model reads 

HLc\i^n>= M^\^n>, (46) 

where Hlc = P^Pn is the invariant mass (squared) operator, refered for convenience 
to as the light front Hamiltonian of positronium and > being the corresponding 
eigenfunction; n labels all the quantum numbers of the state. Projecting this equation 
on the ee state yields 

HlfM{ee)n>=Ml\{ee)^>.. (47) 

In the ee sector, the effective light-front Hamiltonian consists of the free part and the 
effective electron-positron interaction 

Hl'<^ = Hfl + , (48) 
The light front equation Eq. (ETI) is then expressed by the integral equation 



D 2(27t)^ ^ ^^l^^'c'^l^' *2 > ^«(^' *2) = , (49) 



Sl,S2 



Id 2(27r 



where the integration domain D is restricted by the covariant cutoff condition of Brodsky 
and Lepage 

— < + , 50 

x(l — X) 

allowing for states which have a kinetic energy below the cutoff A. The effective ee 
interaction, Vll/ , contains two-point dynamical and instantaneous interactions (with 
the kernel Wee given by Eqs. (^) and (P5|)) as well as the electron self-energy terms 
in ee sector which are renormalized by the electron mass counterterm. We include the 
exchange and annihilation channels in the effective interaction, Vll/ = Vexch + Vann- 

Numerical solution of the integral equation (^91) produces positronium spectrum, 
without (Fig. ^) and with (Fig. |^ annihilation channel. Including the annihilation chan- 
nel, we obtain for the singlet-triplet splitting a supprising agreement with the equal time 
calculations. Degeneracy of triplet state shows that though calculations are done in the 
light-front framework the rotational invariance is manifest. Next excited states {n = 2) 
also show good agreement with the results of the equal time perturbation theory. 
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Figure 6: The invariant mass-squared spectrum Mf for positronium versus the projection 
of the total spin, J^, excluding annihilation with exponential, Gaussian and sharp cutoffs. 
The number of integration points is A^i = A^2 = 21. 
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Figure 7: The invariant mass-squared spectrum Mf for positronium versus the projection 
of the total spin, J^, including annihilation with exponential, Gaussian and sharp cutoffs. 
The number of integration points is A^^i = A^2 = 21. 
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2.1.2 Effective Hamiltonian for the light-front QCD 



QED calculations have been done in terms of bare unrenormalized parameters and success 
of these calculations is based on aplicability of perturbation theory in terms of small QED 
charge. The same naive picture does not hold for QCD; coupling constant is strong. 
However, as we show below (see p| for details), renormalization of the gluon energies 
might change the situation and QCD calculations can be still made using flow equations. 

The effective quark- ant iquark interaction including instantaneous and djTiamically 
generated by flow equations terms is given by 

= -AiTasCf{rY)B^u, (51) 



which formally coincides with the equation (|40| ) for the ee interaction in QED. However, 
the gluon interaction kernel is given by 



Bf,^ = lim 

/xo,Ao— »0 



2 



(52) 



where the integrals Ii and I2 are defined by two similarity form factors containing the 
cut-off dependent momenta transfer in each vertex. 



00 



and with interchange 1 and 2 indeces for 12- The gluon interaction kernel Eq. (|5^)) is 
defined in the limit of zero gluon mass, /iq = 0, and zero renormalization point, Aq = 0, 
that is the renormalization condition in the gluon gap equation (see below). In the light- 
front framework, momenta transfer are given by 

2/,N _ {x'k_i_ — xk'jj'^ + m'^{x — x'y ^ 

Vll^j - ; t^'renK'^ 

XX 

Ql{\) = Ql{\)\x^{l-.x); x'-*{l-x') , (54) 

where firen{^) is the renormalized gluon mass. On mass shell, momenta transfer are 
reduced to an effective gluon energy, 

Q?(A)=Q^(A)-g-^ + /iL(A), (55) 

which contains three gluon momentum transfer and the renormalized gluon mass. Due to 
the presence of the renormalized gluon mass flow equations eliminate off-diagonal matrix 
elements even for vanishing gluon momenta, ^ = 0, that corresponds to degenerate states. 
This cannot be achieved within perturbation theory (PT), since it breaks down at g = 
and an effective interaction has zero energy denominators. Nontrivial solution for the 
effective gluon mass (gluon gap) is a consequence of a special behavior of the light-front 
QCD in the collinear limit, x~^ —>■ 0, in the contrary to QED (there is no mass gap in 
QED). As mentioned before, this behavior of QCD might be attributed to QCD vacuum 
effects. 

The renormalized gluon mass is defined as a solution of the gluon gap equation, /i(A), 
renormalized by the canonical gluon mass counterterm of the perturbation theory, mcT, 

/^L(A) = /^'(A)+m^^(A). (56) 
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One-point flow equations for tlie liglit-front energies provide tlie gap equations for the 
quark (m^) and gluon (/i^) renormalized masses 
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(57) 



where energies are given by 
kl + ii\X) 



QliX) 



x{l — x) 
k^ + w? 
x{l — x) 
kl + m2(A) 



/i2(A) 



+ 



kl 



/i2(A) 



:i-x) 



m\X). 



(5^ 



These equations ( |57| ) are coupled, since unknown cut-off dependent quark and gluon 
masses enter the r.h.s. of both equations. The energy denominators, Eq. (|58D, show a one- 
loop structure. When expanded in the coupling constant at large cut-offs, A = Auy, these 
equations are reduced to the Dyson-Schwinger equation in the rainbow approximation. 
We solve the gluon gap equation with a constant current quark mass, using special 
prescription to regulate coUinear small light-front x in the integral over dx. Absorbing 
trivial perturbative cut-off behavior by a canonical mass counterterm, Eq. (^61), we obtain 
the renormalized cut-off dependent effective gluon mass. 

Taking into account the cut-off dependent renormalized gluon mass, the quark- 
antiquark interaction (Eq. (|^)) has the following form with three choices of similarity 
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(59) 
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that should be compared with the electron-positron interaction in QED, Eq. (^51). On 
the energy shell, the qq interaction is reduced to the sum of the Coloumb, l/(f, and 
linear confining, l/q^, potentials. Here, a plays the role of the string tension depending 
on the cut-off which regulates small light-front x. This result cannot be obtained in the 
perturbation theory by naive summation over perturbative gluon exchanges in a few 
low orders. There are complex calculations, where infinite orders of given subclasses of 
diagramms result in a confining interaction. 
Summary II 

1) Though calculations are performed in the fixed light-front gauge where the light- 
front and perpendicular coordinates are treated differently, we obtain the covariant result 
using flow equtaions (triplet states are degenerate; effective interaction depends on the 
covariant four-momenta and SQ^). Dynamical terms generated by flow quations are 
crucial to maintain covariance. 

2) Flow equations allow to obtain results beyond the perturbation theory. Also, a 
covariant result corresponding to resummation of large number of perturbative diagrams 
may be obtained by flow equations straightforward. 

3) Flow equations allow to track the covariance and gauge invariance much easier than 
in the standard perturbation theory. (To maintain covariance in the PT, all diagrams 
in the given order including crossed diagrams should be summed. Except for the planar 
diagrams summation of high orders of PT is a problem.) 

2.2 Flow equations for QCD in the Coulomb gauge 
2.2.1 Flow equation scaling 

Flow equations eliminate matrix elements which couple states with large energy differ- 
ences, greater than = A, and later more degenerate states. This is reminisent of the 
energy scaling separation underlying perturbative scaling. Flow equations treats different 
energies separately, starting from UV and scaling down towards low energies, separating 
physics at each characteristic scale. 

In strong interactions/QCD, moving from UV to IR, we encounter the following char- 
acteristic scales and corresponding regimes: at Ajjv the conformal canonical perturbation 
theory (PT) without any scale and with bare quarks and gluons as degrees of freedom; 
at Aqcd PT breaks down due to the dimensional transmutation in the renormalization 
group (RG) {Aqcd is the RG invariant scale); at y/a conflnement steps in and hadron 
bound states (instead of quarks and gluons) become relevant degrees of freedom; at 
chiral symetry dynamically breaks, leading to the Goldstone boson-pion; nuclear physics 
decouples further down the scale with and vr as degrees of freedom. Physical ranges 
decouple from each other due to weak interaction (say in the chiral PT, interaction be- 
tween the Goldstone mode and hadrons is weak). Flow equations decouple these regions, 
providing dynamics on all energy scales including crossover between weak and strong 
coupling. 

In what follows, we consider regime where physics is dictated by phenomena of con- 
flnement and chiral symmetry breaking and take into account using flow equations effects 
of high energy QCD which influence the hadron scale physics. 
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2.2.2 Model calculations in QCD 



Calculations made in the Coulomb gauge QCD illustrate clear the scheme of flow equa- 
tions. The question we address is, how starting with canonical QCD Hamiltonian to 
obtain using flow equations an effective block-diagonal Hamiltonian with flxed number 
of Fock components in each block. It seems, that this goal is possible to achieve only for 
a 'conflned QCD' with constituent quarks and gluons. Indeed, in the presence of a strong 
conflning interaction current (bare) quarks and gluons are renormalized, getting dressed 
and becoming constituent quarks and gluons. This QCD motivated Hamiltonian can be 
written in a Fock space of constituent quarks q and gluons g as 



/ qq 



HqcDiq.g) 



V 



qqg 



\ 



99 
qqqq ) 



(60) 



where in every next Fock state a qq pair or a gluon are added. In terms of constituent 
quasiparticles, there is a natural energy gap of order of one GeV between different sec- 
tors. In the gauge fleld theories, mixing between sectors is mediated by minimally coupled 
gauge fields and is strongly suppressed. Small parameter is no longer the coupling con- 
stant, but rather a ratio of the off-diagonal matrix element which mixes sectors to the 
energy gap between diagonal sectors. In the canonical QCD, the off-diagonal sector is 
given by the mixing Coulomb interaction which value is of order of an inverse Bohr ra- 
dius or a current quark mass - several MeV, and the mass gap between sectors is defined 
by a confining scale with a value of order of a constituent mass - one GeV; i.e. a small 
parameter is 
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Ml - M2 



lOMeF 



(61) 



where V\2 is a matrix element of canonical interaction between the first and second Fock 
state with masses Mi and M2 correspondingly. Thus, in the constituent basis, mixing 
between sectors is suppressed and can be eliminated perturbatively by fiow equations 
with a small parameter given by Eq. (0). However, perturbative expansion in this pa- 
rameter holds when working between sectors but not inside a diagonal sector where a 
strong confining interaction dominates. In the effective block-diagonal Hamiltonian, the 
diagonal sectors should be diagonalized numerically. 

At low energies one could think about introducing hadron states as elementary degrees 
of freedom, and constructing a phenomenological Hamiltonian for strong interactions of 
the form 





/ meson 
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Hstrong{hadron) = 




hybrid 






[ 




gluehall J 



(62) 



where hadron states are put on diagonal with increasing glue content from up-left corner 
to down-right; off-diagonal blocks (here empty cells) contain possible interactions which 
mix different diagonal blocks (Fock sectors). In order to find the physical (pure) states 
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one should block-diagonalize this Hamiltonian matrix, ehminating "off-diagonal" inter- 
actions. Probably, one could do it perturbatively, since as known from phenomenology 
the mixing between hadron Fock states is suppressed and a small parameter for iterative 
procedure could be 

M{hybrid) — M{meson) 1 

M{hyhrid) To ^ ^ ' 

where M{hybrid) denotes hybrid mass, etc. This description is along with effective low- 
energy models of nuclear physics, where nucleons are degrees of freedom and pions n 
are mediating interactions between them. 

Summarizing, it seems that the 'confined' QCD provides an effective description 
which matches degrees of freedom in both Hamiltonians given by Egs. and (p^. In 
this framework, using flow equations perturbatively we reduce bound state problem to 
an eigenstate equation with the lowest Fock component of quasiparticles. 



2.2.3 Duality and BV transformation 

Now that we know how to work with the 'confined' QCD which resembles CQM, the 
question is, how to get this effective theory. In other words, how to transform the original 
canonical QCD which is strong coupled, with nonfixed number of particles and has 
complex vacuum with confinement and chiral symmetry breaking to the CQM type 
effective theory which is weak coupled, where the valence quarks and gluons are confined 
in bound states, with the chiral symmetry dynamically broken, and has simple vacuum 
state. We are looking for a dual transformation between these two theories, in the sense 
that the duality is between the strong interacting QCD and weak interacting effective 
theory. 

Many body approach suggests, that the BCS type Bogoliubov-Valatin (BV) transfor- 
mation might fill the gap between QCD and CQM. We adopt that the BV transformed 
QCD has chiral invariant simple vacuum, but instead contains chiral noninvariant in- 
teractions including strong BCS interactions, which reflect confinement and the chiral 
symmetry breaking explicitly, and residual dynamical interactions, which are weak. Ap- 
plying flow equations to the BV transformed 'confined' QCD, we eliminate off-diagonal 
dynamical interactions and obtain the block- diagonal effective Hamiltonian containing 
strong confining chiral noninvariant interactions in addition to dynamical interactions in 
each block. Diagonalizing sector Hamiltonians numerically, we obtain eigenstates which 
are confined and are not invariant under the chiral transformation, and involve dynamics. 
Due to the dynamical interactions, we obtain for the first time equations and correspond- 
ing solutions which are both UV and IR finite. Also, for the first time, solutions include 
explicitly confinement and chiral symmetry breaking as well as dynamically propagating 
gluons. 



2.2.4 QCD motivated Hamiltonian in the Coulomb gauge 

One of the ways to implement confining potential in the QCD Hamiltonian is to use the 
Couomb gauge, which is defined as V ■ A = 0. The Coulomb gauge QCD Hamiltonian 
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includes the free Hamiltonian, Hq, and two types of interactions: instantaneous, W, 
describing static properties, and dynamical, V, involving propagating gluons, 

H = Ho + W + V. (64) 

The free Hamiltonian (quark and gluon kinetic energies) is given by 

Ho = j dxijj^x) {-ia ■ V + Pm) ip{x) 

+ Ti fdx (n\x) + B\{x)) , (65) 



where the non-abelian magnetic field is B = Bi = VjA^. — VkAj + (yf[Aj, A^], and its 
abelian part is represented by B^. Dynamical interaction includes the minimal quark- 
gluon coupling, Vqg, and the non-abelian three- and four-gluon interactions, Vgg, i.e. 
V — Vqg + Vgg. Explicltly they are given by 

= ~9 j dxip"^ {x)a ■ A{x)ip{x) 

Vgg ^ Tr J dx [jU(x)J-^U(x) - W^(x)) +Tr J dx (b\x) - B\{x)) . (66) 

The Coulomb gauge fixing produces the instantaneous quark-quark and gluon-gluon 
interactions, 

W = ^g'ldxdyJ-'p''ix){x,a\iV-D)-\-V'){V-D)-'\y,b)Jp\y), (67) 

with the leading order Coulomb behavior, 1/q'^. Here, the charge density, p, contains both 
quark and gluon components, p"-{x) = ip'^{x)T"-il){x) + /"^'^A^(a;) • n'^(a;). A complete 
solution of the Coulomb gauge constraint equations is encoded in the Fadeev-Popov 
determinant, J = det(V • D) with £> = V — gA, which is unknown. We assume, that a 
nonperturbative solution of J produces the linear confining potential. Thus, we consider 
QCD motivated Hamiltonian, given by the Couomb gauge canonical Hamiltonian with 
the instantaneous interaction including the Coulomb and linear confining potentials; i.e. 
W Wo, 



Wo = ~j dxdyp''{x)VL+ci\x-y\)p''{y) 



a 



CfVL+c{r) = ar-Cf^. (68) 

This QCD motivated Hamiltonian defines our model, used for calculations of glueball 
and meson bound states and presented below. 

2.2.5 Effective Hamiltonian for gluodynamics in the Coulomb gauge 

In a pure ghiodynamics, charge density contains only the gluon component. Following 
the same strategy as before, we represent physical fields, which are three gluon field. A, 
and its conjugate momentum, H, in a second quantized form, 

^"W = /(^7fef^^^'^ + "'"'^-'^^^" 

"^^^^ = (^/yt"'"^^^ " «f (-k)]e^'^- . (69) 
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Since a confining potential has been introduced, a trivial perturbative vacuum does not 
insure minimum for the ground state. Ground state is shifted to some unknown nonper- 
turbative state |0)Arp. The annihilation operator is defined to annihilate in this vacuum, 
i.e. a\0)NP = 0; and creation operator acting on this vacuum produces quasiparticle 
(effective gluon) with unknown energy u!{k). Gluon energy uj{k) is kept as a trial param- 
eter through out the calculations and is found variationally by minimizing the vacuum 
(ground) state energy. The canonical commutation relation is 



[<(k),a7(k')] = (27r)W3)(k-k')A,(k) 



(70) 



where the gluon operators a"(k) = J2\=i,2^i(^^ A)a"(k, A) are transverse, i.e. k ■ a"(k) 
k • a"^(k) = 0; so that polarization sum is 



Ai(k) = -^)ei(k. A) = Sij - kikj 



(71) 



A=l,2 



with the unit vector ki = ki/k; and ki ■ -Dij(k) = 0. 

In a second quantized form, the Coulomb gauge Hamiltonian of gluodynamics con- 
tains the following terms; the gluon kinetic energy is given by 



Hoil) 



+ 



dk 

(2^ 

k2 



+ a;(k,0)af(k)<(k) 



^^^^ -^(k,/))-«(k)<(-k)+h.c.) 



(72) 



The instantaneous gluon-gluon interaction (Eq. (|68|)) is given by 



L+C 



1/2 



yL+c(ki + k2): a,'(ki) + af(-ki 



a-fko 



l)uj{k3,l) 

f(-k3) 



where VL+c(k) is a Fourier transform in the momentum space of linear confining plus 
Coulomb potentials. 



V, 



L+C 



(74) 



here the adjoint Casimir is Cadj = N^. The nonabelian gluon part (Eq. (|66|)) includes in 
the order 0{g) a. triple-gluon coupling 



^3.(0 



L\ (2vr)= 



(2vr)^5(^)(Ek. 



^ijk{ki, k2, k^ 



uj{ki, l)uj{k2, l)uj{k3, 1) 



(75) 



go{ki, k2, ks, /)<(ki)4(k2)a^(k3) + ... 



with 



^ijk{ki,k2,k3,l) 



6 



{{ki - ksjjSik + {k2 - ki)k6ij + {k^ - k2)i6jk) 



(76) 
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In the order 0{g^) the normal-ordered four-gluon vertex reads 



^uj{ki, l)uj{k2, l)uj{k3, l)uj{ki, I) 



(77) 



a,^(ki) +af (-ki) a^(k2) + af(-k2) af (kg) + af (-kg) a'A\i,) + af {-\^^ 



Couphng constants and masses (relevant and marginal operators, respectively) are func- 
tions of the flow parameter / and momenta involved. In addition terms from normal 
ordering with respect to the NP vacuum |0)7vp arise; condensates Oq; Ol+Ci O^g and 
gluon polarization operators Hl+c, ^ig- 

Hamiltonian matrix of gluodynamics is depicted in the Table |^. It includes the triple- 
gluon vertex (for example, element -^14) and the instantaneous gluon-gluon interaction 
(for example, element Hi^); (four-gluon vertex is not depicted). Off-diagonal elements 
mixing different Fock sectors are eliminated by flow equations. 

Elimination of the off-diagonal elements through the second order generates the 
Hamiltonian matrix depicted in the Table ^. Explicit form of the effective Hamilto- 
nian is not given here (for details see 0). Black area in the effective gg interaction (i?|f ) 
depicts the interaction generated by flow equations. 

We solve the effective Hamiltonian in the two lowest sectors. This gives the gap 
equation and the Bethe-Salpeter bound equation for the glueball states. 

Minimizing the ground state energy (Hff), 

»*)=0, (78) 

we obtain the gap equation for the effective gluon energy, uj{k), which is a function of 
gluon momentum k and the cut-off A. Gap equation can be also obtained by demanding 
the off-diagonal one-body operator (i^if), having the structure (aa -|- a^a^), to be equal 
to zero (thus it is not depicted in the Table Explicitly the gap equation is given by 



P + m^y(A) 
dq 1 



+ j'V. / V,,o(k - q) (1 + Ck4f) {^l - 



dq 1 
dq 1 ^(k, q) 



J (27r) ^ t^qt^k-q t^q + U^k-q 

where the mass counterterm mcT(A) is defined by 

m^^(A) = -^N, A\ (80) 

TT O 
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Table 3: Hamiltonian of gluodynamics in the Coulomb gauge. 
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Table 4: Effective Hamiltonian of gluodynamics generated through the second order by 
flow equations in the Coulomb gauge. 
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which absorbs the leading ultraviolet divergences in Hg// as A — oo. In the regulating 
fuction, we assume all external momenta to be soft compared to momentum flowing in 
the loop. Here, the tensor structure of generated terms is given by 

G{k, q) = 2(1 - {kqf) (^k' + q' + + (^9)')) , (81) 

and Vl+c is the Fourier transform of the sum of linear and Coulomb potentials. The 
obtained gluon gap equation is both UV and IR finite. UV divergent behavior arising 
from the Coulomb interaction, four-gluon, and dynamically generated terms is cancelled 
by the canonical gluon mass counterterm. 

We approximate a glueball bound state to consist of two valence constituent gluons. 
The glueball wave function in the rest frame is 

l^n) = / 7TSi</'n^'(q)«f (qX(-q)lo)ivP, (82) 



(2vr) 



where a^(q) creats a quasiparticle from a nontrivial vacuum |0)jvp with an effective 
dispersion relation ci;(q). Since in this constituent basis an effective Hamiltonian He// 
is block-diagonal, the mixing with four gluon and higher states is suppressed. Therefore 
Tamm-Dancoff approach with Hg// is a reasonable approximaion for bound states. To 
obtain the Tamm-Dancoff bound state equation we project the Schrodinger equation 
Heffli^n) = Enl'^An) ou the two-body sector; the result reads 

m[Heff,a^Kq)a^H-q)m = {En-E,) [x^ic^cj^' {c^ , (83) 
with the notation 

(x;/,.,(q)0^^'(q)) = A.'(q)^.y (q)0n^"(q) + d,,,{c^)d,mW{-^) , (84) 

here color indices are omitted, the polarization sum Dij is given in Eg . (|7lD , and -f^e// 
is given in the Table ^ In Eq. ( p3D we subtracted the vacuum energy Eq, defined as 
Heff\0)NP = Eo\0)np- Explicitly Tamm-Dancoff equation for scalar and pseudoscalar 
glueball states with total angular momentum J, parity P, and charge conjugation C and 
only the instantaneous interactions included is given by 




^(p)^(q) 



0n(g) 



V'..c(p - q) '"''''t"'f'' f^""(P.q)4.(p) . (85) 

^(p)^(q) 



with 



F°++(p,q) = 1 + 

i^°"+(p,q) = m)- (86) 
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Here the Coulomb mass counterterm is given by 



Qis ^^2 



(87) 



and Vl+c is the sum of hnear and Coulomb potentials in the momentum space. One-body 
sector of the effective Hamiltonian contributes to the kinetic part of the Tamm-Dancoff 
equation, and two-body effective interaction to the potential part. UV divergences, com- 
ming from the Coulomb potential, are cancelled by the mass counterterm. IR divergent 
behavior, due to the confining potential, cancels in the potential and kinetic parts with 
each other. Such complete cancelation happens only for the color-singlet state, and does 
not occure for colored objects. 

For completeness, we summarize the complete Tamm-Dancoff equation in the scalar 
and pseudoscalar channels when all the terms up to the second order (instantaneous, 
four-gluon and dynamical interactions) are included. It is give by 
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(88) 



with the total angular momentum J, parity P, and charge conjugation C. Here the tensor 



structures are given by 


F°++(p,q) 




i^°-+(p,q) 


= 2(pg) 


i?°++(p,q) 
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p2g,2 



2(p-q)2 



^pq (p - q)2 
Again, this equation is both UV and IR finite. 



(89) 
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q(GeV) 

Figure 8: One particle dispersion relation. Dots represent the numerical solution of gap 
equation c<j(k) (a^ = 0.4, a = 0.18Ge\^^, A — AGeV, Nc — 3), the solid line stays for the 
free dispersion relation a;(k) = k. 
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Figure 9: Gluon mass. Dots represent the numerical solution for m(k) = a;(k) — k {ag — 
0.4, a = O.lSGeV^, A = AGeV, N^^S), the solid line is a fit m(k) = 0.9 * exp(-A;/0.95) 
(parameters are in GeV). 
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Figure 10: Cut-off dependence of the effective gluon mass (a^ 
O.lSGeV^, N^ = 3). 
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Figure 11: Gluon condensate (a^ = 0.4, a = O.lSGeV^, iVc = 3). 

Numerical results 

Numerical solution of the gap equation is shown on Figure ^ with dots for numerics, 
solid line is a perturbative free behavior, u = k. At large momenta an effective gluon 
energy covers the perturbative behavior, while at small momenta u tends to a constant. 
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an effective gluon mass at zero momentum. We obtain m(0) = 0.9 GeV with the string 
tension a = 0.2 GeV^. We define an effective gluon mass as the difference of the effective 
gluon energy and the perturbative value, Figure || (dots is numerics, solid line is a 
parametrization with an exponential fall off). At large momenta, there is zero mass 
and at g ~ we have m(0) ~ 1 GeV. 

Effective gluon mass at zero gluon momentum as a function of the cut-off is given at 
Figure [1^ (in logarithmic scale). Leading behavior is absorbed by the mass countert- 
erm. The logarithmic dependence is left, which is relatively slow. 

Gluon condensate, calculated using the numerically obtained gluon energy, is another 



nonperturbative chracteristic. It grows as a function of A logarithmically (Fig. 11). We 
attribute this slow dependence to a numerical artefact. We regulate the gluon condensate 
by subtracting the perturbative contribution 



TT 



pa jpa \ 
fiu fiu/ 



dkk a. 



(^(k) - kf 



(90) 



Using the dispersion relation obtained above uj(k) the gluon condensate is obtained 
1.3 ■ 10~^GeV^^ (for the cut-off A = 4 GeV), that agrees with the sum rules. (Therefore 
A = A GeV has been used for all our numerical calculations.) 

Numerical calculations of the Tamm-Dancoff equation. Eg . (^5|) , are performed vari- 
ationally with a set of gaussian test functions. Results of calculations for the lowest 
glueball states are presented in the Table (|^) and compared with the available lattice 
data. 





0++ 


0*++ 


0-+ 


0*-+ 


Tamm-Dancoff, (MeV) 


1760 


2697 


2142 


2895 


lattice data, (MeV) 


1730(80) 


2670(130) 


2590(130) 


3640(180) 



Table 5: Glueball spectrum for the lowest scalar and pseudoscalar states (a^ = 0.4, a = 
0.18GeV'^,A = 4GeV,Nc = 3). Lattice data are from C. Morningstar and M. Peardon, 
hep-lat/990100^ . 

Lattice calculations are done for SU{3) pure gluodynamics, using anisotropic lattice 
and improved SII action. Better agreement with the lattice data is achieved for the scalar 
channel. Remarkable, the mass of the lowest scalar glueball O"*""*" is roughly twice of the 
effective gluon mass m(0) obtained before. This confirms the constituent picture. The 
better agreement is achieved for the scalar masses, indicating that dynamical terms are 
important for the excited states and the complete Tamm-Dancoff equation with all terms 
should be solved. 



2.2.6 Effective Hamiltonian for QCD in the Coulomb gauge 

This model is also suitable to study the chiral symmetry breaking. Now, quarks are 
dynamical and the charge density has only the quark component. The instantaneous 
interaction contains the Coulomb plus linear confining potentials. Gluons are pertur- 
bative and quarks have unknown dispersion relation, dynamical energy E{k) which is 
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parametrized through the BV angle $(fc). Vacuum \Q) contains condensates of quark 
pairs (i.e. it is the BCS vacuum). The Fock space is constructed from this vacuum using 
quasiparticle operators and which appear in the field expansions 

= E / (^[^(^' ^)^(^' ^) + ^(-^' ^)^^(-^' '^y^"" 

with b\Q) = d\Q) = 0. Here the quark operators are given in the helicity basis and all 
descrete numbers (helicity, color, and flavor for the quarks and color for the gluons) are 
collectively denoted by s and a, respectively. We use spinors in the massive basis, where 
nonzero effective quark mass explicitly is included in the spinor 



( 



1 



u{k,s) = V^(fc) + M(fc)(^^.^/(^(^)^^,(^)JX. 



l + sik)xs 




- s{k){(T ■ k)xs 



1 - s(fc)(cr ■ k){-icr2Xs 



1 + s{k){-ia2Xs 
with sine and cosine of the Bogoliubov angle $(fc) given by 

sin($(fc)) = s{k) = , ^^^^ , cos(<l>(fc) = c{k) - ^ 




(92) 



k' + AP{k) Jk'^ + M^{k) 



E{k) = ^Jk' + M^{k), (93) 

E{k) is a single-quark energy which we refer as a gap energy below. Effective quark mass 
M{k) is kept as an unknown variational parameter through out the calculations, and is 
found from the gap equation by minimizing the ground state (vacuum) energy. Gluon 
energy is uj{k) = k. 

With the definition Eq. (p2D, the spinors satisfy the nonrelativistic normalization and 
orthogonality relations u'^ {k, s)u{k, s) = v^{—k,s)v{—k,s) = 1 and u''{k, s)v{—k, s) = 
v'^{—k, s)u{k, s) = 0. Canonical (anti) commutation relations are 

{b{k, s), b\k', s')} = {d{-k, s), d\-k', s')} = {27if 6{k - k')6s,s' 

[a,{k, a), a]{k', a')] = {27rf 6{k - k')D,j{k)6a,a' , (94) 
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where the gluon operators a = aiik)"- = X]a=i,2 ^ j(^5 ^^e transverse, i.e. k ■ 
a°'{k) = k ■ a°'^{k) = 0, and Dij{k) is a polarization sum 

A=l,2 

with unit vector component ki — ki/\k\ and ki ■ Dij{k) — 0. 

Following the same scheme, we represent the Hamiltonian in a second quantized form. 
The free Hamiltonian (kinetic quark and gluon energies) is given by 

+ {ks{k, I) - mc{k, l))[bl{k)dl{-k) + d,(-k)6,(fc)]) 

+ J 7^s^ik)at\k)am ■ 



(27r) 

The quark-gluon vertex is given by 



(96) 



^o(fei, k2,ks, l)ds, {-ki)T%,, (fe2) ^i&L (97) 

V2^(fc3) 

vl(-ki)aiu,,{k2){27rf5^'\ki -k^- k^) + ...1 , (98) 



where E{k) = yk"^ + M'^{k) and Lo{k) = k. By implementing flow equations effective 
coupling constants are generated {go{ki, k2, k^; I) ...) which are functions of all three 
momenta corresponding to a given Fock sector and depend upon the flow parameter I. 
The instantaneous interaction includes the linear confining and Coulomb potentials and 
is given in the quark sector by 

HL+c{i) = E / f n ^] (27r)^5(^)(fci + ks-k2- k,)VL+c{k^M) (99) 

S1-S4"' \n=l V^^J / 

: [<(fci)6t^(fcO + <(-fei)rf.i(-fci)]r»K(fe2)&..(fc2) + ^..(-fe2)<(-fe2)] 

K(fe3)&l3(fe3) +<(-fc3)43(-fc3)]T»K(fc4)6.,(fc4) +^.4(-fe4)4(-fc4)]: , 

where VL+c(fc, O) ^ VL+c{k - q), 

CfVL^cik) = 27rC^^ + 47r^, (100) 

with the fundamental Casimir operator Cf = T"'T°' = {N^ — l)/2Nc = 4/3. The spinors 
are also functions of the flow parameter, i.e. everywhere Us{k, I) and Vs{k, I). In addition, 
the condensate Oq, Ol+c and self-energy terms Ei+c arise due to the normal-ordering 
of the canonical QCD Hamiltonian with respect to the BCS vacuum. Dynamical quark 
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\qq > 


\qqg > 


\qqqq > 




|0 > 














\q > 














jq> 












\qqg > 












\qqqq > 













Table 6: QCD Hamiltonian in the Coulomb gauge. 

energy is expressed through the BV angle, which is kept as a trial parameter through 
calculations. 

Matrix elements of the Coulomb gauge normal-ordered Hamiltonian are shown in 
the Table For example the quark-gluon coupling is given by the element ifi4 and the 
instantaneous interaction is given by H^^. Off-diagonal elements are eliminated by flow 
equations. 

As a result of this elimination new terms arise in the diagonal sector, depicted in the 
Table 0. 

The gap equation allows determination of a nontrivial vacuum with quark condensates 
and propagating quasiparticles, here quarks with a dynamical mass. There are several 
ways to obtain this equation, the most common based upon a variational principle to 
minimize the vacuum (ground state) energy. The variational parameter is the angle of 
transformation from undressed to dressed particle (quasiparticle) operators, which 
defines a quasiparticle basis (Eqs. (|9T|) and (|92|) ) with a dynamical quark mass M{k). 
Therefore, analogosly to the gluon sector, minimizing the vacuum energy of the effective 
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qqqq > 
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IQ > 
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\qqg > 












\qqqq > 
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Tabic 7: Effective QCD Hamiltonian generated through the second order by flow equa- 
tions in the Coulomb gauge. 
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Hamiltonian, 

we obtain the gap equation for the unknown $(fc) or M{k). Using the condensate terms 
given by Hi(^ , the gap equation is obtained as 



ks{k) - m(A)c(fe) = / JS^^CfVL^cik, q) [c{k)s{q) - s(k)c(q)k ■ q 



(27r) 
dq 
(27) 

where I — k — q, the potential functions are given by 



+ / j^,CfW{k,q)[c{k)s{q)-s{k)c{q)k-iq-i 



CfVL+c{k,q) = o7ir^ + 



2(fe-q)2 (k-q)^ 

CfW{k,q) = -J-, 71 7T, (103) 

^ ^ u;(k-q)(E{q)+u;(k-q))' ^ ' 

and the running mass m(A) includes the mass counterterm and is defined by 
m(A) = m + Mct{^) =m(l - -^^6 In A ) 




Mcr(A) = -6m = -j^6mlnA. (104) 



{Any 

The mass counterterm is proportional to the bare quark mass, m, and thus vanishes in 
the chiral limit m ^ 0. Quark mass counterterm cancels UV divergencies. In the chiral 
limit, no counterterms are required. The gap equation in UV finite due to the dynamical 
interactions generated by flow equations. 

Next we consider the quark condensate for a single quark flavor. We regulate 

the quark condensate by subtracting the perturbative contribution 

mm-W'pm = -A'c/7^(Tvsi"(*)-Tvsf w) 

where m is the bare quark mass. As |fe| oo the mass gap M{k) m and the non- 
perturbative sine behaves as s{k) — >• m/\/W+ln?. Thus this subtraction improves the 
convergence of the quark condensate integral in the UV. Here the equal-time propagator 
IS given by S^^\k) = 1/20 [mA -j-k{l + B)]. When the scalar part of the propaga- 
tor is nonzero {mA ^ 0), the nonzero mass gap (Mq = mAjiX -|- S)) and the chiral 
condensate ((■(/'■?/') o ~ / dkmA/Q) are generated. 

We approximate the meson bound state as a valence pair of quark and antiquark, that 
corresponds to the Tamm-Dancoff approximation. In terms of the quasiparticle operators 
the TDA meson creation operator reads 

Ri = j ^,j:bl{q)d\{-q)^l^i\q) . (106) 
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Figure 12: The numerical solution of dynamical quark mass, Mo{k), and the scalar part 
of the propagator, 2k(7s, in the chiral limit when Coulomb and generated potentials with 
the running coupling as{k'^) are added to confinement. The parameters for the numerical 
solution of the gap equation are a — O.lSGeV^, A = IGeV. The results are compared 
with the fit function given by Mo{k) = 0.0060/(A;2[ln(A;V0.04)]0-^3) 



acting on the vacuum, it creates a meson with a wavefunction \ipn) with quantum number 
n; annihilation operator defines the nonperturbative vacuum, i.e. 



Rn\^) 



0. 



(107) 



Projecting Schrodinger equation i/e//|'*/'n) = En\ipn) onto one particle-one hole truncated 
Fock sector, we get the TDA equation 

{Q\[R^,[H,ff,bi{k)dl{-k)m = M^^fik), (108) 

where the binding energy is defined as M„ = — Eq. Commutation relation of the 
meson operators 



{n\[Rn>,Ri]\9)^N5nn', 
leads to a normalization condition for the wave functions 



I 



(109) 



(110) 



with the normalization constant A'^. 

In the Random Phase Approximation, we allow in addition to the qq creation the qq 
annihilation in vacuum, with corresponding wave functions X and Y for each correlation. 
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Figure 13: One particle dispersion relation, Eo{k) = \Jk'^ + Mo(k)'^, free dispersion, 
Eo{k) = k, and the vector part of the propagator, l/2cr^, in the chiral limit with con- 
finement and with confinement plus perturbative potentials (i.e. when Coulomb and 
generated potentials are added). The parameters are same as in Fig. |T2. 



Corresponding generalization of the operator given by Eq. ( |106|) containes meson creation 
and annihilation terms 



The RPA wavefunction and the RPA vacuum are defined by 

Qn\^) = 0, (112) 

where, though the same notations were used above, they should not be confused with 
the TDA wave function and TDA vacuum. The RPA bound state equation is given as a 
system of equations for both components of the wave function, X and Y, 

{Qm,,[H,ff,bUk)dl{-k)]]\Q) = M^Xfik) 

{Q\[Q^,[H,ff,b^{k)dp{-k)]]\Q) = M^Yf{k), (113) 
From the meson commutation relation 

iSl\[Q^,,Qi]\n) = N5^n', (114) 

the following normalization condition for the wave function components, X and Y , is 
obtained 



dq 

(2^ 



E [x'j\q)K\q) - Yi?*{q)Y^"'{q) ] = N6^^, , (115) 
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Figure 14: Sine and cosine of the Bogoliubov-Valatin angle in the chiral hmit with con- 
finement and with confinement and with confinement plus perturbative potentials (i.e. 
when Coulomb and generated potentials are added). The parameters are same as in Fig. 

0- 

with the normalization constant A^. 

We considered psedoscalar tt- meson (J^*" = O"*""*" L = S = 1, J = 0) and vector p- 
meson {J^'^ = 0~~^ L = S = J = 0) channels to illustrate effects of the chiral symmetry 
breaking. Based on quantum numbers of these states, the tensor structure of vr- and 
p-wave functions can be identified as 

Xf{k) = {za2r^X^{k),Y:^{k) = i-za2r^Y^{k) 

Xf{k) = {aia2r^X,{k),Y;^{k) = {-ia2(Tr^Yp{k). (116) 
Normalization is chosen 

^ {X*{k)Xik) - Y*{k)Y{k)) = 1 , (117) 
reducing the normalization costant of the full wave function to 

iV= (^„|^n) =2iVe, (118) 

where factor 2 comes from trace in the spinor space, and Nc = 3. The RPA equations 
for the momentum wave function components X{k), Y{k) have the same form for vr and 
p states 

M^Xik) = 2e{k)X{k) - I ^^4.(fc, q)X{q) - J ^^4,(fc, q)Y{q) 
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q^dqdx 
47r2 



GUk,q)X{q) 



q^dqdx 
47r2 



-M„r(fc) = 2e{k)Y{k)-j 
(fdqdx 



q^dqdx 



G,y{k,q)Y{q) 

' q^dqdx 



{k,q)Y{q)-J^-^Iy^{k,q)X{q) 



I 



47^2 



Gyy{k,q)Yiq) 



q^dqdx 
47r2 



(119) 



where the kernels / and G for tt are given by 



7:,(fe,qr) = 7;rfc,q) = C/T4+c(fc,q)-[(l + s(fe))(l + s(q)) 



+ (l-s(fe))(l-s(9)) + 2c(fe)c(q)a:l 

1 
2 

- (l-s(fe))(l + s(g)) + 2c(fe)c(g)x: 



= j;,(fc,q)=C;Vl+c(fc,q)-[-(l + Kfe))(l-^(q)) 



= Gl(k,q)^2CfW,{k,q)- 



(l-s(k))(l + s(q))-2c(k)c(q) 

1 



{l + s{k)){l-s{q)) 

(l + x'^)kq - x{P + g2) 



= G;,{k,q)^2CfW,{k,q)^ 



+ {l-s{k)){l-s{q))-2c{k)c{q) 
and for p are given by 



{k-qy 
{1 + s{k)){l + s{q)) 

(1 + x2)/cg — x{k'^ + q^) 



{k - qf 



(120) 



KAk.q) 



I'Jk,q) = CfVL+c{k,qh 



{l + s{k)){l + siq)) 



1 



+ -(1 - s{k)){l - s{q)){Ax^ - 1) + 2c{k)c{q)x 



P(k,q) = PMq)^CfVL+cik,q)- 



2 

- hl-sik))il + siq)) + lc{k)ciq)x 



-i(l + 5(fe))(l-.(9)) 



GS,(fc,q) = G';^(k,q) = CfW,{k,q) 



,{1 + s{k)){l - s{q)) 



+ \{1 - s{k)){l + .(q)) (l - ^) - |c(fe)c(q) (x + 



(1 — x'^)kq' 



GL{k,q) = G'^,(fc,q) = C;iy2(fe,q)- 



2 L3 



(fe - g)^ 
^{1 + sik)){l + siq)) 



ik-q) 



+ ^(1 - s{km - s{q)){2x' - 1) + ^c(fe)c(g) [x - 



(121) 



where x — k ■ q] and we have used 

(l-a;2)g2 



l-(fc-Z)^ 



(fe-qr) 



2 ,l-(9-0 



2 (l-a;2)fc' 



(fc - qr)^ 
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TDA, (MeV) 


RPA, (MeV) 


conf. 


504 


1364 


2115 


222 


1416 


2298 


conf.+Coul. 


608 


1514 


2249 


427 


1521 


2309 


conf.+Coul.+gen. 


513 


1411 


2161 


180 


1413 


2218 



Table 8: Pion spectrum for the ground, first and second exited states in the TDA and 
RPA approaches with confining, confining+ Coulomb and confining+Coulomb+generated 
potentials taken. Chiral limit m = (a^ = 0.4, a = 0.18^6^^, A = lOGeV^). 



k ■ Iq ■ I 

with I = k — q. 



x{P + q^) — (1 + x^)kq 



:i22) 




Figure 15: Cut-off dependence of the constituent quark mass in the chiral limit (same 
parameters as in Fig. |12]). Crosses represent solution with confinement. Boxes [diamonds] 
represent solution when Coulomb and generated potentials are added with the constant 
value of coupling Qs [with the running coupling gs{q^)]- 



Numerical results 

For the numerical calculations we have used the routines from the SLATEC linear 
algebra archive, part of the Netlib database maintained by UTK and ORNL. The routines 
are found at www.netlib.org/slatec/lin/ and a description of the entire SLATEC archive 
can be found at www.netlih.org/slatec/toc. 

Effective quark mass and energy are shown in Figures |12| and 0, respectively. At 
large momenta effective mass vanishes and the energy tends to the perturbative value, 
E{k) = k, while for zero mometum, the energy tends to the constant, which we refer as 
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Figure 16: Quark condensate cut-off dependence in the chiral limit (same parameters 
as in Fig. |T^. Crosses represent solution with confinement. Boxes [diamonds] represent 
solution when Coulomb and generated potentials are added with the constant value of 
coupling Qs [with the running coupling Qsiq"^)]- 



m, (MeV) 


TDA, (MeV) 


RPA, (MeV) 


150 


1038 


1926 


2936 


1037 


1986 


3077 


100 


885 


1762 


2697 


868 


1811 


2826 


50 


716 


1590 


2431 


660 


1626 


2537 


10 


553 


1446 


2212 


366 


1460 


2283 


5 


532 


1428 


2186 


293 


1437 


2250 





513 


1411 


2161 


180 


1413 


2218 



Table 9: Pion spectrum for the ground, first and second exited states in the TDA and RPA 
approaches for different current masses of constituents. Confining+Coulomb+generated 
potentials are taken, the same parameters as in the Table |. 



the constituent quark mass. No counterterms are required in the chiral limit. Due to the 
dynamical interactions we recover the large momentum behavior correctly, correspond- 
ing to the UV asymptotic region, which behavior is obtained in the course of explicit 
numerical solutions of the QCD renormaliztion group and operator product expansion 
and QCD sum rules. 

At low momenta cosine behaves lineraly (Fig. 0), with the slope 1/M(0) (inverse of 
quark constituent mass), c(fc) = k/M{0). 

The sensitivity of the constituent quark mass, M(0), and the chiral condensate in 
the chiral limit to the cut-off A is displayed in Figures |l^ and 16. Only when all terms. 



44 





TDA, (MeV) 


RPA, (MeV) 


conf. 


659 


1484 


2258 


642 


1482 


2256 


conf.+Coul. 


750 


1678 


2515 


732 


1676 


2514 


conf.+Coul.+gen. 


718 


1592 


2377 


700 


1590 


2376 



Table 10: Spectrum of the p meson for the ground, first and second exited 
states in the TDA and RPA approaches with confining, confining+ Coulomb and 
confining+Coulomb+generated potentials taken. Chiral limit m = 0, the same parame- 
ters as in the Table |[ 



m, (MeV) 


TDA, (MeV) 


RPA, (MeV) 


150 


1130 


2086 


3247 


1128 


2086 


3247 


100 


986 


1916 


2990 


983 


1915 


2990 


50 


839 


1744 


2692 


833 


1744 


2692 


10 


727 


1616 


2436 


714 


1615 


2435 


5 


719 


1603 


2406 


704 


1601 


2405 





718 


1592 


2377 


700 


1590 


2376 



Table 11: Spectrum of the p meson for the ground, first and second exited states 
in the TDA and RPA approaches for different current masses of constituents. 
Confining+Coulomb+generated potentials are taken, the same parameters as in the 
Table |. 



instantaneous and dynamical interactions are added, we obtain the stable result. This 
proves that the gap equation is renormalized completely. We need to add the renormal- 
ization group running of the coupling to renormalize the condensate. For the first time, 
completely renormalized stable mass and condensate are obtained in the chiral limit, 
with values Mo(0) = IQMeV and {'4>i))o = -(155 MeV')^. Flow equations improve the 
chiral condensate by 68%. 

Relative scales of constituent masses are obtained: for the gluon < IGeV, for the 
s-quark > 200 MeV, for the u, rf-quarks < 100 MeV. 

Results for the pion and p-meson masses in RPA and TDA approaches are presented in 
Tables |8|-pT|. In the chiral limit RPA gives ground state pion masses which are significantly 
lower than those obtained using TDA. Including the dynamical interaction terms reduces 
ground state pion mass even more, encreasing the mass splitting between the pion and 
p meson (Tables |] and |10|). In the chiral limit we get 



M^ = 180 MeV , Mp = 700 MeV 

Mp - M^ = 520 MeV . (123) 

The IT — p mass splitting of 520 MeV in the chiral limit is close enough to the lattice 
data splitting of 600 MeV. However, we are unable to get zero mass pion either in the 
BCS or adding the leading order fiow equations corrections. One of the reasons is break 
of the covariance in this model. Flow equations improve the vr — p mass splitting by 32% 
in the TDA and by 24% in the RPA. 
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Figure 17: 71 — p mass splitting for the ground state in TDA and RPA approaches. 
Chiral hmit, m = 0, and confining+Coulomb+generated interactions are taken (same 
parameters as in the Table |^). Fit function is 12000/(m+20). 

1.09 -1 1 



M(A)/M(10) 
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Figure 18: Cut-off dependence of the pion mass in the chiral limit M(A), normalized to 
the pion mass at A = lOGeV, M(10), (same parameters as in the Table H). Line with 
crosses represents RPA solution with confinement + Coulomb, the line with boxes stays 
for RPA solution when generated potentials are added. 
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Numerically obtained dependence of the n—p mass splitting, Mp—M^, as a function of 
the bare mass of one of the quarks is shown in Figure \L7\. We find l/rriconst behavior, where 
TTT'const IS a coustitueut quark mass, which is valid for heavy quarks and continues to be 
valid for lighter constituent quarks. This fall of is more rapid in the RPA than in the TDA 
approach. From the RPA fit function, a constituent quark mass can be approximated, 
uniformly for heavy and light quarks, as rriconst = m+20 (MeV), with the bare quark mass 
m. The 1/ rriconst behavior is characteristic for the hyperfine interaction. We find in the 
chiral limit, that in the TDA roughly 30% of the vr— p mass splitting is due to the presence 
of the hyperfine interaction and 70% due to the chiral symmetry breaking. In the RPA 
this ratio is 40% for the hyperfine and 60% for the chiral symmetry breaking. However, 
the numerical value of this ratio depends on the details of the confining interaction. 

The dependence of the pion mass in the chiral limit on the cut-off parameter, M(A), 
is shown in Fig. 0. The RPA solution, with only confining+ Coulomb potential included, 
grows unlimited (due to Coulomb), while adding the generated term stabilizes M(A), 
which saturates roughly at M(10). Stable result confirms that the TDA/RPA equations 
are completely renormalized when the generated by fiow equations terms are included. 
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Figure 19: Ground state pion wave function in TDA and RPA approaches. Chiral limit, 
m = 0, and confining+Coulomb+generated interactions are taken (a^ = 0.4, a = 
0.18Ge\/^A = lOGeV). 



Pion wave functions for the ground and first excited states are depicted at Figures 19 



and pO. Using the TDA and RPA wave functions, we obtained the pion decay constants. 



Though we obtained low value of quark condensate —(155 Me V^)^, we obtain a realistic 
pion decay constant, 92MeV, in the RPA. 
Summary III 

1) For the first time, due to the dynamical interactions, the UV and IR finite equations 
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Figure 20: Pion wave function for the first exited state in TDA and RPA approaches. 
Chiral hmit, m = 0, and confining-|-Coulomb-|-generated interactions are taken into 
account (same parameters as in Fig. |l^. 

are obtained in many-body technique. In the chiral hmit, no additional UV renormaliza- 
tion is required (no momentum dependent counterterms) . 

2) For the first time, the instantaneous as well as dynamical interactions are included 
(In this way we have utilized the success of the BCS model, where nonperturbative 
features such as dynamical chiral symmetry breaking and massive quasiparticle modes 
are explicitly present, and included perturbatively dynamical interactions in the BCS 
framework). In the RPA roughly 40% of the vr — p mass splitting is due to the presence 
of the hyperfine interaction and 60% due to the chiral symmetry breaking (link between 
different model calcualtions) . 

3) Due to the dynamical interactions, closed Ward-Takashi identities are obtained. 
Though we work in the gauge theory, there are no ambiguities, since all fields in the 
effective Hamiltonian are physical. 



3 Conclusions 

In this work we have outlined a strategy to derive an effective renormalized Hamiltonian 
by means of flow equations. Application of the flow equations with the condition, that 
particle number conserving terms are considered diagonal and those changing the particle 
number off-diagonal led as in other cases to useful effective Hamiltonians. 

The main advantage of flow equations as compared to other many-body approaches 
is, that states of different particle number are completely decoupled, since the particle 
number violating contributions are eliminated down to A = 0. Thus one is able to truncate 



48 



the Fock space and the bound state problem reduces to an eigenstate equation in the 
lowest Fock sector of quasiparticles. Positronium (meson) problem is approximated as a 
state of two valence electrons (constituent quark and antiquark), described by a Bethe- 
Salpeter equation with an effective Hamiltonian in two-body sector. Energy/mass of 
a quasiparticle is found from the gap equation with an effective Hamiltonian in one- 
body sector. Coupled gap equation and Bethe-Salpeter equation with effective sector 
Hamilonians are analyzed further analytically and numerically for a one-particle mass 
gap and bound state spectrum, respectively. 

We applied flow equations to the light-front Hamiltonians of QED and QCD and 
to the Hamiltonians of gluodynamics and QCD (with dynamical quarks added) in the 
Coulomb gauge. We found that dynamical new terms generated by flow equations are 
crucial in order to obtain an agreement with covariant calculations and the experimental 
results. In the light-front positronium problem, due to dynamical terms triplet states 
are degenerate and thus rotational invariance (that means covariance in this case) is 
maintained. Using flow equations, we obtain the singlet-triplet splitting in the light-front 
framework which agrees with the experimental data. 

In the Coulomb gauge QCD, due to the dynamical interactions we obtain for the first 
time the gap equation and the Bethe-Salpeter equation which are finite both in the UV 
and IR regions. We find that no additional renormalization is required in the chiral limit. 
For the first time, we include dynamical interactions in the framework of many-body 
technique where confinement and chiral symmetry breaking explicitly present. In this way 
we have studied an overlap of hadron physics and high energy QCD. We took into account 
the hyperfine interaction as well as chiral symmetry breaking and obtained the tt — p 
mass splitting caused by the instantaneous and dynamical interactions. We find in the 
RPA, that roughly 40% of the n — p mass sphtting is due to the presence of the hyperfine 
interaction and 60% due to the chiral symmetry breaking. Flow equations improve the 
gluon and quark condensates as compared with other Hamiltonian methods (we obtain 
(a/7rF^t.F^i.) = 1.3 ■ 10~^GeV agrees with the sum rules and {4>4>) = —{l55MeV)^ is 
still low). 

We find that for glueball states confinement plays the dominant role and dynami- 
cal terms are necessary to insure the correct (finite) UV behavior. However, in meson 
systems, dynamical terms are equally important with the instantaneous terms responci- 
ble for chiral symmetry breaking. Using fiow equations, we obtain masses of low lying 
glueballs and mesons which agree with lattice data and the experimental results. 
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The main publications on this topic are: 
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